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A diffuse flux of astrophysical neutrinos above 100 TeV has been observed at the IceCube Neutrino 
Observatory. Here we extend this analysis to probe the astrophysical flux down to 35 TeV and 
analyze its flavor composition by classifying events as showers or tracks. Taking advantage of lower 
atmospheric backgrounds for shower-like events, we obtain a shower-biased sample containing 129 
showers and 8 tracks collected in three years from 2010 to 2013. We demonstrate consistency with 
the (/ e : fy, : / T )© « (1 : 1 : 1)® flavor ratio at Earth commonly expected from the averaged 
oscillations of neutrinos produced by pion decay in distant astrophysical sources. Limits are placed 
on non-standard flavor compositions that cannot be produced by averaged neutrino oscillations but 
could arise in exotic physics scenarios. A maximally track-like composition of (0 : 1 : 0)© is excluded 
at 3.3cr, and a purely shower-like composition of (1 : 0 : 0)© is excluded at 2.3<r. 


Introduction —Traveling virtually unimpeded through 
matter, radiation, and magnetic fields, astrophysical neu¬ 
trinos probe otherwise inaccessible regions of the high- 
energy universe. If produced from cosmic rays interact¬ 
ing with gas and radiation at their sources, they convey 


unique information about astrophysical particle acceler¬ 
ators in their direction, energy, and flavor [1-5]. Though 
no individual sources of TeV cosmic neutrinos have yet 
been found, a diffuse flux was observed by the IceCube 
Neutrino Observatory above 100 TeV in three years of 
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data [6, 7]. Here this work is expanded to observe the 
diffuse astrophysical neutrino flux down to 35 TeV and 
derive constraints on its flavor composition. 

Astrophysical neutrinos are expected from the decay 
of secondary particles such as pions, kaons, muons, and 
neutrons produced in cosmic ray interactions. In the 
model of diffusive shock acceleration, the differential en¬ 
ergy spectrum of injected cosmic rays follows a power 
law oc E ~ 7 with 7 ~ 2 [8, 9]. Though this spectrum may 
be modified by propagation in cosmic magnetic fields en 
route to Earth [10], neutrinos produced at the source 
retain the same spectral index 7 as the injected cosmic 
ray spectrum, provided the environment is sparse enough 
to allow particles to decay rather than interact. In the 
most commonly considered scenario, the decay of pions 
and their daughter muons dominate the neutrino flux, 
resulting in a flavor ratio of (/ e : : f T )s = (1:2: 0)s 

at source [11, 12]. However, the composition could vary 
from (0 : 1 : 0)g to (1 : 0 : 0)5 under a multitude of sce¬ 
narios including muon energy loss in high matter density 
or magnetic fields [13-16], muon acceleration [17], and 
neutron decay [18]. 

As first noted in [19], neutrino oscillations, averaged 
by propagation over astronomical distances, transform 
the flavor ratio according to the PMNS mixing matrix 
[20-22]. Taking global best-fit mixing parameters [23], a 
flavor ratio at Earth of ( f e : f^ : f T )@ = (0.93 : 1.05 : 
1.02)0 « (1 : 1 : 1)0 is expected for a (1 : 2 : 0)5 source 
composition, a result of the near tribimaximal form of the 
PMNS matrix [24], For a composition at sources varying 
from (0 : 1 : 0)s to (1 : 0 : 0)s, the composition at Earth 
varies linearly from (0.6 : 1.3 : 1.1)0 to (1.6 : 0.6 : 0.8)0. 
Though expected to be negligible [11], even a large u T 
contribution at sources causes only a small deviation from 
this range. Because of this limited variation in the fla¬ 
vor ratio at Earth for all possible source compositions, 
the observation of a ratio inconsistent with these expec¬ 
tations would signal new physics in the neutrino sector, 
such as neutrino decay [25, 26], sterile neutrinos [27], 
pseudo-Dirac neutrinos [28, 29], Lorentz or CPT viola¬ 
tion [30], and quantum gravity-induced decoherence [31]. 
Measuring the flavor ratio of astrophysical neutrinos is 
interesting both as a probe of the source of high energy 
cosmic rays and a test of fundamental particle physics. 

Neutrino interactions in IceCube —The IceCube detec¬ 
tor consists of 5,160 digital optical modules (DOMs) in¬ 
strumenting lkm 3 of clear ice at the South Pole [32, 33]. 
Each DOM contains a photomultiplier and digitizing 
electronics that detect Cherenkov light emitted from sec¬ 
ondary particles produced in neutrino interactions [34]. 
Neutrino events are generally classified into two topolo¬ 
gies: track-like, where the path of an outgoing charged 
particle is visible, and shower-like, where the region of 
light emission is too small to be resolved and appears 
point-like. For both topologies the energy deposited 
within the detector can be reconstructed within ~ 15% 


above 10 TeV [35]. Neutrino direction can be recon¬ 
structed with a median angular error of < 1° for tracks 
versus ~ 15° for showers above 100 TeV [35]. 

In charged-current (CC) interactions, a neutrino de¬ 
posits its energy into a charged lepton and a hadronic 
shower, and the topology of an event depends on flavor. 
For v e CC interactions, the outgoing electron initiates 
an electromagnetic shower indistinguishable from the ac¬ 
companying hadronic shower. For z/ M CC interactions, 
the outgoing muon leaves a long track in addition to a 
hadronic shower. If the muon leaves the detector, the 
deposited energy is only a lower bound on the neutrino 
energy. For v T CC interactions, the outgoing tau decays 
quickly and is difficult to resolve for energy < 1 PeV 
[19, 35, 36]. However, tracks may be observed from the 
muonic decay of the tau with 17.4% branching ratio. In 
neutral-current (NC) interactions of all flavors, a neu¬ 
trino deposits on average ~ 1/3 of its energy into a 
hadronic shower but with a cross section ~ 1/3 of the 
CC cross section [37]. Above ~ 10 TeV, neutrino fluxes 
become attenuated by interactions in the Earth, though 
v T fluxes are regenerated by subsequent tau decay to neu¬ 
trinos [38]. 

The backgrounds in astrophysical neutrino searches are 
muons and neutrinos produced by cosmic ray air show¬ 
ers in Earth’s atmosphere. Muons dominate the trigger 
rate in IceCube and usually create long tracks. How¬ 
ever, they can also appear shower-like if they undergo a 
single catastrophic energy loss inside the detector. At¬ 
mospheric neutrinos are usually divided into two cate¬ 
gories. The first arises from the decay of kaons, pions 
and muons, producing mostly v jL [39-41]. Since time 
dilation causes decay to be less likely than interaction 
at high energy, the neutrino energy spectrum is asymp¬ 
totically one power steeper than the primary cosmic-ray 
spectrum, and the angular distribution is peaked at the 
horizon. Time dilation also suppresses v e from muon de¬ 
cay down to the detector depth, and the remaining v & are 
from kaon decays at the level of k, 4%. The flux of 
atmospheric v & has recently been measured in the TeV 
range by IceCube [42], 

The second category, yet to be observed, results from 
the prompt decay of charm mesons [43-48], yielding a 
nearly equal mixture of v e and but negligibly small v T 
[49]. Henceforth referred to as the charm neutrino flux, 
it should follow the same ~ EV 2 7 spectrum as primary 
cosmic rays and also be isotropic. Since atmospheric 
backgrounds from muons and produced in ir/K de¬ 
cay are largely track-like, astrophysical events dominate 
over backgrounds down to lower energies in the shower 
channel, and a contribution from charm decay may be 
more easily identified [50]. 

Event selection- Data collected at IceCube in 974 days 
from May 2010 to May 2013 are used in this analysis. 
During the design of the event selection criteria, 90% of 
data was kept blind. Following the same strategy as the 
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FIG. 1. The log likelihood ratio between shower and track 
reconstructions for veto-passing events with more than 1500 
photoelectrons. Error bars are 68% Feldman-Cousins inter¬ 
vals. The contribution of muons is determined from a muon 
control sample. The dotted lines show the total amount of 
Vn CC events (pink) and all non-z^ CC events (maroon) from 
the best-fit distributions of astrophysical and n/K neutrinos. 
The last bin contains all overflow events with LLR > 500. 


previous 3-year analysis [6, 7], events are selected using 
an outer layer of DOMs to veto the vast majority of in¬ 
coming muons, isolating neutrino interactions of all types 
starting within the detector from across the entire sky. 
Also similarly, the muon background rate is estimated 
with a control sample. Using outer DOMs to tag incom¬ 
ing muons, an inner volume geometrically similar to the 
full fiducial volume is defined with its own veto layer of 
DOMs. After correcting for its smaller fiducial volume, 
the rate of tagged but unvetoed events yields the muon 
background rate in the full detector. 

Down-going atmospheric neutrinos can also be vetoed 
by accompanying muons from the same cosmic ray air 
shower that reach the detector [51]. The veto probability 
is determined using the analytic calculation described in 
[52], accounting for muons from both the same decay as 
the neutrino and other decays in the air shower. The re¬ 
sulting suppression of down-going atmospheric neutrino 
events distinguishes them from the isotropic distribution 
expected from a diffuse astrophysical flux. For the charm 
neutrino flux, otherwise isotropic, this suppression is the 
only distinguishing feature if the astrophysical flux has 
a power-law index close to 2.7 and a non-standard flavor 
ratio (1:1: 0)©. 

Showers and tracks are classified by performing a per- 
event maximum likelihood analysis of the first photon 
arrival times in every DOM. Each event is reconstructed 
according to the hypothesis of an infinite track with 
constant light emission along its path [53] and the hy¬ 
pothesis of a point-like shower [54], yielding likelihoods 
L '©rack and Tshower- A log likelihood ratio, LLR = 


— ln(Lghower/LTrack) 5 is formed, with negative values be¬ 
ing considered showers and positive values tracks. 

Figure 1 shows a distribution of LLR for veto-passing 
events in the 10% unblind data sample producing more 
than 1500 total photoelectrons (PE). Best-fit neutrino 
distributions for a (1 : 1 : 1)© composition (discussed 
later) are shown, and the prediction for muons is ob¬ 
tained from the corresponding likelihood ratio distribu¬ 
tion in the muon control sample. The agreement with 
data illustrates that the control sample reliably predicts 
the rates of both shower-like and track-like muons. Also 
shown is the combined distribution of astrophysical and 
atmospheric CC events, illustrating that most are clas¬ 
sified as tracks. The remaining 30% of CC events clas¬ 
sified as showers arise when the outgoing muon has too 
little energy to be resolved or escapes near the edge of 
the detector. 

In the final selection, shower-like events above 1500 
PE are selected, while for tracks, only events above 6000 
PE are selected due to the larger background from pen¬ 
etrating muons. The deposited energy and direction of 
each event is reconstructed using the full timing distri¬ 
bution of recorded photoelectrons in every DOM. For 
events classified as showers, point-like light emission is 
assumed, whereas for tracks, the energy deposition is un¬ 
folded along the path of the track [35]. To avoid system¬ 
atic uncertainties relating to muons, down-going, shower- 
like events below 20 TeV are excluded because they are 
dominated by muons. 

Statistical Analysis —To measure the flavor ratio of the 
astrophysical flux, we follow the approach of earlier Ice- 
Cube analyses [55] and perform a binned, maximum like¬ 
lihood fit over the 2D distributions of deposited energy 
and reconstructed declination of both showers and tracks. 
The expected count in each bin is calculated from Monte 
Carlo simulation and depends on a set of nuisance pa¬ 
rameters, which describe systematic uncertainties, and 
physics parameters, which are of interest to be measured. 
The likelihood contains two terms—the first describing 
the Poisson distribution of bin counts and the second 
penalizing deviations of nuisance parameters from their 
central values according to their uncertainty [56]. To con¬ 
struct confidence intervals and perform hypothesis tests, 
we use the profile likelihood method [57] and minimize 
the likelihood over nuisance parameters. Unless noted, 
we assume that profile likelihood ratios follow a y 2 dis¬ 
tribution [58]. 

Systematic uncertainties are either detector-related, 
such as DOM optical efficiency and ice optical properties, 
or theoretical, such as neutrino fluxes and cross sections. 
All tend to uniformly scale the rates of tracks and show¬ 
ers, maintaining their ratio and leaving their energy and 
angular distributions unaffected. Thus, nuisance param¬ 
eters describing backgrounds become overall rate scaling 
factors applied to the distributions of i t/K neutrinos, 
charm neutrinos, and muons. To describe atmospheric 
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FIG. 2. The best-fit deposited energy distributions for showers and tracks, divided into the southern (down-going) and northern 
(up-going) samples, assuming a power-law astrophysical flux with (1:1: 1)© flavor ratio at Earth. Showers in the southern sky 
below 20 TeV are dominated by muons and excluded; however, the prediction from the control sample measurement is shown 
in this region. Though not shown, 4 bins in declination ( 8 ) are used in the fit. 


neutrinos, we use the flux calculation of HKKMS [59] 
for n/K neutrinos and ERS [60] for charm neutrinos. A 
correction is included to describe the cosmic ray knee in 
the model of [61], and the HKKMS flux is extrapolated 
above 10 TeV as described in [55]. No priors are used 
to constrain the scalings of these atmospheric neutrino 
distributions. For muons, although the control sample 
constrains the overall passing rate, there are insufficient 
statistics in its energy and angular distribution, so sim¬ 
ulation is used. It is based on a parametrization of the 
deep-ice muon flux [62, 63] obtained from CORSIKA air- 
shower simulation [64] and the cosmic ray model of [61]. 
Since muons are the dominant background for astrophys¬ 
ical tracks, we allow the scalings applied to track-like 
muons and shower-like muons to float independently, ac¬ 
counting for uncertainties such as ice properties, energy 
loss cross sections, and muon bundle multiplicity that 
could skew the ratio of tracks to showers. These scalings 
are, however, constrained by a Gaussian prior of 8.4±4.2 
events each, derived from the 4 surviving events each in 
the track-like and shower-like muon control samples. 

When placing limits on the flavor ratio, nuisance pa¬ 
rameters also include those describing the astrophysical 
flux. For this analysis, we use an isotropic, power-law 
flux with the following parametrization for each neutrino 
flavor, 

$«(£) - 3$ 0 /a,® ( 10 o TeV) ’ ^ 

where / a ® is the fraction of each flavor at Earth, 7 is 
the spectral index, and $0 is the average flux of v and 
v per flavor at 100 TeV. An equal v and v flux is as¬ 


sumed. Though this does not hold for neutrinos of photo- 
hadronic origin, there are consequences for a flavor ratio 
measurement only from yet-unobserved v e interactions at 
the 6.3 PeV Glashow resonance [65, 66], too high in en¬ 
ergy to have a significant impact with currently available 
statistics. 

Results- After all selection criteria, 129 showers and 
8 tracks remain in the final event sample, forming a su¬ 
perset of the earlier 3-year sample with 28 showers and 
8 tracks [7]. Before attempting to constrain the astro- 
physical flavor ratio, it is necessary to verify that the 
adopted isotropic, power-law model of the astrophysical 
flux adequately describes the data. Assuming a flavor 
composition of (1 : 1 : 1)® at Earth, the best-fit distri¬ 
butions are shown in Fig. 2. Noteworthy are the best- 
fit astrophysical flux parameters, 7 = 2.6 ± 0.15 and 
$0 = (2.3 ±0.4) x 10 -18 GeV -1 s -1 cm -2 sr -1 . While 
being compatible with the previous 3-year result [7], the 
spectral index is substantially different from 7 = 2, which 
is rejected at 3.0er. The preference for 7 > 2 comes 
mostly from low-energy data rather than a lack of events 
above several PeV. A high-energy cutoff in the astro- 
physical flux of the form a E~ 2 exp (—E/E c ) is also dis¬ 
favored with respect to a power law at 2.9<r. Finally, the 
power-law model with (1:1: 1)® flavor ratio is in agree¬ 
ment with data with a goodness-of-fit p-value of 0.13. 

A large charm flux to explain low-energy data is dis¬ 
favored since the suppression of down-going events by 
accompanying muons is not observed, and the best-fit 
scaling of the ERS flux is 0. Even fixing the ERS scaling 
at its 90% upper limit of 3.4 obtained here, the astro- 
physical index only changes to 7 = 2.5. These results 
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FIG. 3. The exclusion regions for astrophysical flavor ratios 
(/e : /a* : /t)® a t Earth. The labels for each flavor refer 
to the correspondingly tilted lines of the triangle. Averaged 
neutrino oscillations map the flavor ratio at sources to points 
within the extremely narrow blue triangle. The « (1 : 1 : 1)® 
composition at Earth, resulting from a (1 : 2 : 0)s source 
composition, is marked with a blue circle. The compositions 
at Earth resulting from source compositions of (0 : 1 : 0)s and 
(1 : 0 : 0)s are marked with a red triangle and green square, 
respectively. Though the best-fit composition at Earth (black 
cross) is (0 : 0.2 : 0.8)®, the limits are consistent with all 
compositions possible under averaged oscillations. 


are consistent with a recent, dedicated IceCube measure¬ 
ment of the astrophysical spectral index and charm flux 
with improved veto techniques [62]. Nuisance parameters 
describing n/K neutrinos and muons are also consistent 
with expectations from the HKKMS flux and the control 
sample measurement. 

This analysis is sensitive to the astrophysical flux in the 
neutrino energy range 35TeV - 1.9 PeV. The lower and 
upper bounds of this range, E \ ow and E up , were calcu¬ 
lated separately by fixing the astrophysical spectral index 
and normalization at their best-fit values, excluding the 
flux with E < £) ow or E > E up , respectively, refitting the 
data with nuisance parameters left free, and finding the 
values for £) ow or E up that decreased the log likelihood 
by 1/2 each. 

With a power-law astrophysical flux describing the 
data, we then further allow the flavor composition at 
Earth to float and calculate exclusion regions according 
to the Feldman and Cousins approach [67], as shown in 
Fig. 3. Though the best-fit composition is (0 : 0.2 : 0.8)® 
at Earth, the limits are compatible with all standard fla¬ 
vor compositions possible under averaged neutrino oscil¬ 
lations at < 68% confidence level. 

With showers and tracks serving as the only two identi¬ 
fiers for three flavors in this analysis, there is an inherent 


degeneracy in the determination of astrophysical flavor 
ratios. This is reflected in the strong anti-correlation 
between v e and v T components, which both produce 
mostly showers. The degeneracy is broken mainly by two 
effects—the shift in the v T deposited energy distribution 
caused by invisible energy lost to neutrinos in tau decay 
and the lack of observed v e Glashow resonance events. 
The preference for a j/ r -like signature is not statistically 
significant, and future work to identify v T signatures at 
PeV energies may resolve this degeneracy. 

Since compositions produced by averaged neutrino os¬ 
cillations (narrow blue triangle in Fig. 3) are nearly or¬ 
thogonal to the flavor degeneracy in IceCube, constraints 
on source flavor composition are possible but not yet sig¬ 
nificant. After restricting to flavor ratios allowed by aver¬ 
aged neutrino oscillations, no source composition can be 
excluded at > 68% confidence level, and this remains true 
even with the additional constraint f T g = 0 expected at 
astrophysical sources. 

Having found agreement with the predictions of aver¬ 
aged neutrino oscillations, constraints are placed on non¬ 
standard flavor compositions producing a large v e or 
fraction at Earth. A maximally track-like, pure sig¬ 
nature of (0 : 1 : 0)® is excluded at 3.3cr and a purely 
shower-like v e signature of (1 : 0 : 0)® at 2.3cr. 

These results contrast with an earlier analysis of Ice- 
Cube’s 3-year data, which found a preference for (1:0: 
0)® over (1:1: 1)® at 92% confidence level [68]. We at¬ 
tribute this discrepancy mainly to two unaccounted for 
effects - partial classification of CC events as show¬ 
ers and systematic uncertainty on muon background. Re¬ 
peating their analysis but accounting for the ~ 30% of v p 
CC events classified as showers and using a profile like¬ 
lihood incorporating the 50% uncertainty in muon back¬ 
ground, a (1 : 0 : 0)® best-fit is still obtained but neither 
(1:1: 1)® or our best-fit of (0 : 0.2 : 0.8)® are ex¬ 
cluded at > 68% confidence level. Since only shower and 
track counts were analyzed, the tighter constraints re¬ 
ported here result from the use of energy and directional 
information in addition to the lower energy data. 

Future measurements of the flavor ratio at IceCube will 
use improved veto techniques, include up-going tracks 
starting outside the detector, and search for high-energy 
signatures of v T . With these improvements, measuring 
the flavor composition at astrophysical sources and pre¬ 
cision tests of neutrino oscillations over astronomical dis¬ 
tances will be in reach. 
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